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Zusammenfassung 

Es wird gezeigt, dass man die Losung gewisser linearer stochastis- 
cher Differentialgleichungen in Hilbertraumen, und zwar soldier mit 
beschrankten Operatoren sowie der in [10] und [15] betrachteten kon- 
servativen stochastischen Schrddingergleichung, analog zur Lie-Trotter- 
Produktformel aus der Halbgruppentheorie gewinnen kann, indem man 
die Gleichungen in einen "deterministischen" und einen "stochastischen" 
Anted aufspaltet und die zugehorigen Ldsungsfliisse abwechselnd iterativ 
auf die gewiinschte Anfangsbedingung anwendet. 
Betreuer: Prof. Dr. Detlef Diirr 



Abstract 

We prove that the solution of certain linear stochastic differential equations 
in Hilbert spaces, namely those with bounded operators as well as the 
conservative stochastic Schrodinger equations analyzed in [10] and [15], 
can be obtained - along the lines of the Lie- Trotter product formula from 
semigroup theory - by splitting the equation into a "deterministic" and a 
"stochastic" part and alternately applying the corresponding solution flows 
in an iterative manner to the initial value. 
Advisor: Prof. Dr. Detlef Diirr 
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CHAPTER 1 



Introduction 



The main subject of this work is the stochastic Schrodinger equation 

(1.1) dip t = {-iH - \A 2 )^ t dt + A^tdtt. 

for a Hilbert space valued diffusion process ipt with standard Wiener process & and, in general, 
unbounded self-adjoint operators H and A. This is an extension of the usual Schrodinger 
equation 

(1.2) dipt = -iHiptdt 

by a stochastic term that will make the wavefunction "collapse" . Note that solving the isolated 
"stochastic part" 

(1.3) dVt = Afadtt ~ \A 2 ^ t dt 

of this equation is considerably simpler than solving the whole equation - in an abstract 
framework as well as if one wants to find an explicit solution formula. The same is true for 
the Hamiltonian part (1.2). Therefore I wondered whether one can find an alternative access 
to the solution of the whole equation by constructing it out of these two partial solutions. 
Such a method would reflect the intuitive notion leading to equations like (1.1) that the two 
mechanisms (1.2) and (1.3), when combined additively, act simultaneously on tpf 

Let me remind you that the analogous question for deterministic equations of the form 

(1.4) x{t) = (A + B)x{t) 

is settled by the Lie- Trotter product formula, which, in the presence of unbounded operators, 
holds true e.g. in the following form (see [3], [16]): 

Theorem 1.1. Let A, B and C := A + B be generators of contraction semigroups, noted 
as e tA etc., on a Banach space. Then s — \im n ^ >00 (e™ A e™ B ) n = e tc holds true with "s-lim" 
denoting the strong operator limit. 

Roughly speaking, this means that one gets the solution of (1.4) by subjecting the initial value 
alternately for infinitesimal time intervals to the solution flows e tA and e tB which belong to 
the single operators. Chernoff has proven the Trotter formula in a more general framework, 
replacing F(t) := e^ A e^ B by an arbitrary contraction valued mapping that fulfills F(t) = 
id + tC + o(t) for t -> 0: 

Theorem 1.2. Let F be a strongly continuous mapping from [0, oo) to the set of contraction 
operators on a Banach space such that F'(0), the clousre of its strong derivative at 0, generates 
a contraction semigroup. Then s — lim^oo F (^) n = e tF '^ . 
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Introduction 



Besides the classical Trotter formula, this theorem comprises e.g. the cases F(t) = (id — tC) 1 
and F(t) = [(id - tA)(id - tB)]- 1 . 

Our goal is to prove similar results for linear stochastic differential equations (SDE). In 
chapter 2, after having taken a look at the solution of the isolated stochastic part of (1.1), 
we develop the theory of weak solutions of linear dissipative SDE - a class of equations 
comprising (1.1) - going back to [10] and mention the results of Mora und Rebolledo ([15]) 
about their regularity. Chapter 3 begins with a generalization of Chernoff 's theorem to linear 
SDE which have matrices or bounded operators as coefficients. A result of this kind does not 
seem to have been published so far. Then we combine our proof with ideas form [10] and 
[15] concerning the handling of unbounded operators and arrive at a product formula for the 
stochastic Schrodinger equation (1.1). 

Shortly before we finished this work, Smolyanov et al. have republished the article [8], where 
they also find a product formula for stochastic Schrodinger equations. However, they admit 
only Hamiltonians arising as the Weyl quantization of a classical phase space function and 
bounded collapse operaors instead of our unbounded operator A. 

In chapter 4 we apply our product formula to two prototypes of quantum mechanical collapse 
models, namely the time-continuous QMUPL model ([5]) and the discrete GRW model ([6]), 
and shows that the first can be constructed as a suitable continuum limit of the latter. 



CHAPTER 2 



Solution theory for linear SDE with unbounded operators 



1. SDE with one operator 

Let £i be a Brownian motion on a probability space (CI, P). A simple example for an infinite 
dimensional SDE is the equation 

, . d/ t (x) =xft(x)d£ t 

[ ' ' fo(x) =g(x). 



for a process taking functions on M as values. Reading this equation as a famliy of ordinary 
SDE indexed by x, each one of which describes a real- valued process ft(x), one obtains the 
solution ft(x) = e x ^~2 x t g{x) (a family of geometric Brownian motions, indexed by x as 
diffusion constant). With regard to more complicated examples, the question arises whether 
such processes belong to certain function spaces, which would make them accessible methods 

of functional analysis. As for our example, e x ^~2 x2t = e~ 2 ( x ~^t) 2+ ^t is bounded as a function 
of x, so g 6 L P (R) implies f t 6 L P (M) a.s., and also smoothness properties of g are preserved. 
However, in the Hilbert space case p = 2 we will be concerned with in the following, the 
computation 

E\\f t \\ 2 = -±=je-£j f 2 (x)dxdy = -J= J J e -kiv-**)*+** g\ x )dydx = j e tx2 g\x)da 

shows that ft remains in L 2 (Q, L 2 (M)) only for very particular initial functions g (e.g. such 
with compact support). Otherwise, also the continuous dependence on the initial data cannot 
be expressed in terms of E|| • || 2 . One can get rid of the troubles with the term e tx by 
considering the process ft(x) := e~^ 1+c ^ x f t (x) with c > 0. This modified process solves the 
equation 

(2 2) d/ t (x) =xf t (x)di t - (1 + c)x 2 f t (x)dt 

fo(x) =g{x). 

The preceding observation suggests us to consider the equation 
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with c > rather than dipt = Aiptd^t also if we are dealing with more general self-adjoint 
operaors A on a Hilbert space H (the restriction to T>(A 2 ) will become clear). The formally 
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expected solution 
(2.4) 



ipt = e 



£ t A-(l+c)tA< 



is well-defined according to the functional calculus because f(t,x,y) := e v x -( l + c ) tx2 ; as a 
function of x with fixed values t > 0, y G R or i = y = 0, is bounded on the spectrum of A. ip t 
lies in L 2 (Q, H) and depends continuously on an arbitrary given initial value ip: For a proof 
let fi^ be the corresponding spectral measure. Then one computes 



E\\ip t 



I 2 =E 



J f 2 (t,x,tt)w{dx) = -jL= j e » J f 2 (t,x,y)fi i ,(dx)dy 



1 



v^vrt 

< y i^(dx) = /^(r) 

- note in particular that in the case c = even 

(2.5) E||rf = |H| 2 



holds - and 



(2.6) 



E||Vt-VII 2 = 



1 



2vrf 
1 + e 



_ 2L 

e 2t 



- 2e 



(l+e)fcc 2 



)//,/, (da;) 



(since the integrand is majorized by a constant). By similar calculations on makes sure that, 
for t > 0, f(t,A,y)ip is differentiable once w.r.t. t and twice w.r.t. y with the expected 
results. For example, 

2 



h 



(f(t + h, A, y)^ - f(t, A, y)$) - [-(1 + c)A 2 f(t, A, y)\ 



J p (/(* + h > x > y) ~ /(*> x > f ) + ^ 2 (! + c)x 2 f(t, x, y)) 2 ^(dx) 

e -(l+c)/ix 2 _ J 



= 2j/x-2(l+c)ta; 2 



+ {l + c)x 2 I ^(dx) ^>0, 



because the estimate |e cx — 1\ = c J* e ct dt < c dt = cx valid for c, x > via 

/ -{l+c)hx 2 - \ 2 



2 / e 



2j/x-2(l+c)tx 



<e 2 J /x-2(l+c)te 2 2 



/l 



+ (1 + c)x 2 



( e -(l+c)hx 2 _ ^2 



+ 2(1 + C) 2 X 4 < 4(1 + c )2 x 4 e 2^-2(l+ C )te 2 



leads to a bounded and therefore ^^,-integrable majorant. 

According to the Ito formula, (2.4) satisfies, at least for t > und e G (0,i), the equation 

ik = il> e + J Atfj s d^ s - J + c^^ds 



2.2 General linear SDE with constant coefficients 
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and it remains to be checked whether the terms on the right converge to their counterparts 
in (2.3) for e — > 0. ip e — )• ifi has been shown in (2.6) and for if) £ T>(A 2 ) we have, since in this 
case e £tA-{i+c)tA 2 CO mmutes with A and A 2 , 

n f \l + c)A 2 ^ s ds - f {\ + c)A 2 ^ s ds\\ 2 <e ^ E\\(l + c)A 2 i; s \\ 2 ds 

JO 2 Je * Jo 2 

= (- + c) 2 e [ e E\\e^ A -( 1+c > A2 A 2 4,\\ 2 d S <(- + c) 2 e f E||^|| 2 ds -> 
2 Jo 2 j 

and 

E|| f Aip s d£ s - [ A^ s d£ s \\ 2 = I Ep^ s || 2 ds< I E||^|| 2 ds^0. 

JO Je Jo Jo 

Let me emphasize once more that the conditions c > and ip £ T>(A 2 ) which are required in 
the abstract ansatz can turn out to be far from optimal in concrete special cases - equation 
(2.2) is trivially solvable for arbirary c and even non-measurable initial values if one reads it, 
as performed for (2.1), as a family of ODE. 



2. General linear SDE with constant coefficients 



Let ^, be independent Wiener processes on (f2, P), ($t)t>o the standard extension 

of the filtration generated by them, J = cr(dt,t > 0), fj a separable C-Hilbert space and 
K,Li,L2, ■ ■ ■ operators on ft with dense common domain T>. We now develop the solution 
theory for the equation 

oo 

d^ t =Y,LjAd{;i-Kip t dt 
(2-7) j=1 

V>o =i> e s), 

for an fj-valued process ipt{w) auf Q with particular focus on the special case of the stochastic 
Schrodinger equation 

(2.8) dip t = {-iH - ^A 2 )^ t dt + Afrtet. 

Even in the case of bounded operators, no explicit solution formula for such equations seems 
to exist, but then at least the existence and uniqueness is ensured by global Lipschitz and 
growth conditions. In the general case, an approximation by bounded operators, as performed 
in [10], at least leads to a weak solution. We present this construction and continue with the 
results of [15] on the regularity of these solutions. 

In the following we require the "dissipativity condition" 

(2.9) ^2 \\Ljipf - 2Re (Kip, ijj) < c||^|| 2 fiir alle ip £ V und ein c > 

j 
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for the operators. It has the role of an a-priori growth estimate: Since a formal application 
of Ito's formula yields 

de- ci ||^|| 2 = 2e- c *^Re fa, L^dtf + e~ ct (-c||^|| 2 - 2Re fa,K^t) + ll^tf) dt, 

3 

V v ' 

<0 

it ensures that, if a sufficiently smooth solution ipt of (2.7) exists, e _rf ||V>t|| 2 is a supermaringale 
and, in particular, E ]]?/;* || 2 < e c *||-0|| 2 gilt. Regularity is only analyzed for the special case 
c = 0. Such equations, obviously comprising (2.3) from the last paragraph as well as the 
stochastic Schrodinger equation (2.8), are called "conservative". 

2.1. Weak solutions. In the sequel, by "weak solution" of (2.7) we mean a process ipt 
satisfying 

(</>^t) = ^,^) + J2 I (L*^^s)de s - [ (K*<t>,i> s )ds 
■ Jo Jo 

P-a.s. for all in a dense set T>* (yet to be specified) on which L*- and K* are defined. As 
for the existence of such a set, note that K* are in fact densely defined because K c := K + | 
is accretive according to (2.9) and thus closable, K = K c — | and K* D K. We assume the 
existence of a core T>* for K* such that all L* are defined on it and 

oo 

(2.10) :J2\\ L *j^\\ 2 <°°- 

J'=l 

Theorem 2.1. In i/ie situation described above there exists a weak solution of (2.7) with 

(2.11) EHVtH 2 < HVIIV*. 

If K + ^id is maximal accretive, the solution is unique under this condition. 

Basic facts about accretive operators can be found in [11] or [17] and will be used without 
proof. Moreover, we need the following possibility of restricting the choice of test functions 
when dealing with weak convergence in L 2 (Q): 

Lemma 2.1. Let be independent Wiener processes on (Q, 5", IP), (3t)t>o the standard 

extension of the filtration generated by them and J = a($t,t > 0). We define recursively 
according to 




stochastic processes, aj denoting arbitrary indicator functions of the form X[a ■,&•]• Then the 
random variables 1$ := 1, 1/ := I;(oo) (i.e. the integral in the recursive definition runs over 
the whole of M + ) form a total subset of L 2 (Q,$); the random variables Ii(t) form a total 
subset of L 2 (f2,3t). 

Proof. According to a version of multidimensional Wiener chaos decomposition de- 
scribed in [18], all iterated integrals of the form J °° f^- 1 ■ ■ ■ f* 1 . . . , t m )&Ct\ • • • d^™ 
with <3? £ L 2 (R m ) form a total subset; since iterated integrals are continuous in the integrand, 



2.2 General linear SDE with constant coefficients 
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the choice of <& can be restricted to the total subset of L 2 (M m ) constituted by products of 
the form X[ai,6i](*i) ' ' ' X[a m ,6 m ](*m)- This proves the first part of the statement. In view of 
Ii(t) = E(Ii | 3t) and the L 2 -continuity of the projection E(- | fo), the Ii(t) form a total subset 
ofL 2 (tt,&). □ 



Beweis von Satz 2.1. For technical simplicity we confine ourselves to m G N instead 
of countably many Wiener processes. We start by building modified Yosida approximations 
of Lj and K by bounded operators: For this sake, let K* be an m-accretive extension of 
K* (K* = K* + § is accretive because K c is) and R n := (id + ^K*)' 1 (n G N). Then 
R n is contractive, s — lim n _ 5 . 00 R n = id and s — rim n _>oo R* n = id (because, applying the 
general formula B* A* D (AB)* to the mutually inverse operators R n and id + ^K*, one 
obtains i?* = (id + ^if c ) _1 , X c denoting an m-accretive extension of K c ). L" := LjR n and 
i^T n := R* n KR n are bounded (they are closed and defined on the whole of Sj), unlike the usual 
Yosida approximation they satisfy the dissipativity condition (J2j ||-^"V'I| 2 — 2Re (K n ip,ijj) = 
J2j \\Lj(Rnip)\\ 2 ~ 2Re {K(R n 4>),R n i>} < c\\R n ^f < c\\ip\\ 2 ) and for all <fi G V* it holds 

true that L™*0 " - ~ > °°> L*0 (because, for </> G V* , U-*4> = RnL*j4> and s - limi?* = id) and 
.f\" n *</> — )■ Kef) (because K n *(f) = R* l K*R n (f>, K*R n isis the usual Yosida approximation of K* , 
so it converges to K* on V*, then i?* — >■ id and \\Rn\\ < 1 imply the statement). As a 
consequence, the equations 



(2.12) 



have a unique solution satisfying 

(2.13) E||V>ni 2 < e ci 



This estimate shows that the ipf are uniformly bounded as mappings from a time interval 
[0, T] to L 2 (Si,Sj). We would like to show equicontinuity and apply the Arzela-Ascoli theorem. 
Since we have no clue concerning the limiting behaviour of expressions like E||L"Vf || 2 > we pass 
over to weak convergence in L 2 (Q,S)) = L 2 (fi, C) <g>fj, i.e. we choose £ G L 2 (f2,C) and <\> G T>* 
and consider the complex-valued functions f n (t) := E£(0, ^"). Due to the boundedness of 
ip?,we only have to consider £ and from dense subsets of the unit speres, therefore we 
restrict to G T>*; since L 2 ($7,C) and Sj are separable (the first statement follows from the 
Wiener chaos decomposition), even an appropriate sequence )zeN suffices. The f n (t) are 
uniformly bounded according to the Cauchy-Schwarz inequality and (2.13) and equicontinuous 
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according to 

m rt 



p t pt 

\r\t) - r( S )\ 2 <Eiei 2 -ek^c -c)i 2 = nJ2 / <*,w>d#- / <<^ n c>dT| 

<2E|J^ /"*(£"* 0,C)d^| 2 + 2E| f {K n (j),^)dT\ 2 
<2^ fn{Lf^)\ 2 dT + 2\t-s\ fmK^^fdT 

j Js Js 

<2|t-s|^||Lf c/)\\ 2 sup E||Vv|| 2 + 2|t- s| 2 ||i^>|| 2 sup E| 

j re[s,t] re[s,t] 

<2\t -s\J2 \\L*H 2 e cT M 2 + 2\t - s\ 2 \\K n 0\\ 2 e cT \\ij T \\ 2 



3 

<CU T \t - s\" -« 2 



Therefore, there exists a subsequence (n\)k e ^ of (n) n6 ^ such that E^ 1 (^i,^" fe ) converges 

uniformly in t, a subsequence (n\) of (ra£) such that also E£ 2 (</>2, ^™ fc ) converges and so 

on. For the diagonal sequence (n|) =: (n&) obviously E£j(<£j, converges for all j € N 
uniformly in t; since 

L 2 (0,ji)9?j^ lim E(i),^)eK 
is a continuous functional, there exists a process ^ with 

lim E^,^)=E^,Vt) 

for all £,0 uniformly in t. (2.13) and the weak closedness of balls imply that this process 
satisfies (2.11). 

In order to show that this is the wanted solution, the convergence of the remaining summands 
in the "weak form" of (2.12), i.e. in 

m pt pt 

E^,vr fc )=w,^)+E^,v / Lf^de s )-m<p, K n ^d S ) r esp . 

~[h J Jo 



(2.14) E£<0, = E£(^> + f\ L T^ n s k )^l - E tf/ K 



t 



(strong Bochner integrals and their stochastic analogue allow for dragging bounded operators, 
in particular scalar products, inside), to their counterparts in (2.7) has to be checked. To 



2.2 General linear SDE with constant coefficients 
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begin with, 

|E£ I {K nk *<p^ n s k )ds-El f (K*<f>, <0 s )ds| = | f El{K nk *d>, ^ k )ds - [ E£(iT>, Y> s )ds| 
Jo Jo Jo Jo 

<| f El(K nk *<t)- K*4>,^ k )As\ + \ [ El(K*(/),^ k - V s )ds| 
Jo Jo 

<(E|e| 2 )2||^ nfc 0-^>|| / (E||^ s rifc || 2 )2ds + ...^0; 

Jo 

for the summands in the middle of (2.14) we restrict £ WLOG to the iterated stochastic 
integrals Ii(t) specified in lemma 2.1. This restriction to a total subset is justified by the fact 
that, according to 

E| ^V?*^" fc >d^| 2 = j\{L n k *^: k )\ 2 ds < \\q k *ct>f j\Us k \\ 2 ds < CTe cT U f, 

the sequences of stochastic integrals which are checked for convergence are bounded. Then, 
similar as above and in addition using the Ito formula, the independence of the Wiener 
processes and the fact that, according to their recursive definition, all Ii(t) are L 2 -martingales 
and have mean value 0, we conclude 



\EI l+1 (t) f(L; k >,e)^-I/, + i(t) !\l*^ s )ah\ 
Jo Jo 



= |E f a^I^Lf^^ds -E f a j {s)I l {s){L*ct>^ s )ds\ 
Jo Jo 

<| l\~^)W){L nk *(t> - L*^^)ds\ + | ^ Eaj(s)Ii(s)(L*(f),tpg k - ii s }ds\ 
Jo Jo 

<(E|/Ks)| 2 )5 \\Lf > - L*^.|| J* E||VC fe || 2 d S + ...->(). 

Now let ipj and ipf be two weak solutions of (2.7) satisfying the same initial condition and 
(2.11) and define 5t ■= ipl — tpf. For the unique solvability of (2.7) it is sufficient to show 
E//(t)(0, 5t) = for all 4> &T>* and l\ already specified. Since So = 0, 

(2.15) ~W){<\>M) = (<f>,S t ) = V f (L*<t>,5 s )dfii - f{K*<t>,8 s )ds 

j Jo Jo 

and, according to the product formula for Ito differentials, 
(2-16) 

W)(^S t )= [ I W)^ L *^ 6 ^- I W)(K*0,5 s )ds 

Jo j Jo j Jo 

+ / Va^s)(L*0,<5 s >ds 
Jo 
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for I > 1. For the first time in this work, the integrands of the arising stochastic integrals 
contain products of two L 2 -processes, so they are in L 1 according to Cauchy-Schwarz, but in 
general not in 1? . Therefore, the integrals are no longer defined in the usual 1? framework via 
the Ito isometry, but only in the generalized sense described e.g. in [9], for which continuous 
paths are already sufficient. In general, such integrals are only local martingales. According 
to Holevo, the uniform boundedness of the Z^-norms of the integrands imply their martingale 
property. I have not found out why this is so - if it is true, then, taking mean values in (2.15) 
and (2.16) and defining mi(t) := KIi(t)5t, one gets 



(2.17) (4>,m (t)) = - I (K*<f>,mo(s))d S 

and for I > 1 



o 



(2.18) (0,m,(t)>= f (Vo 3 ( S )L^,m H ( S ))d S - f {K*<f>, m,(s))d*. 

Jo j Jo 

Now let us prove inductively mi = 0, which will of course imply Eij (t) ((f), St) = 0. Starting 
with / = 0, we somehow have to bring the two terms mo in equation (2.17) together. For this 
sake, we take the Laplace transform, perform partial integration on the right-hand side und 
thus get rid of the additional integral in front of one m®: 

t ft ps 1 ft 



e" As (0,m o (s))ds = - / e~ Xs (K*<f>, m (r))drds = -- / e~ As (K*(P, m (s))ds 
o Jo Jo A Jo 

Multiplying with A and dragging out the scalar products, we get 

(2.19) {\<j> + K*<f>, [ e- Xs m (s)ds) = 0, 

Jo 

(the integrals still being understood as Bochner integrals). Now let A > |. According to the 
assumptions, K + | and thus also K* + | = K + | is m-accretive, so Im (K* + A) = fj and 
therefore (K* + \)D* is dense in S). In view of that, (2.19) implies 

t 

e~ Xs m (s)ds = 



o 

for A > |. The Laplace transform is holomorphic, so this holds for all A > 0; since a function is 
uniquely determined by its Laplace transform, mo = follows. In the induction step l — l—tl, 
the first integral in (2.18) is according to the induction hypothesis, so mi = follows from 
this equation along the same lines as mo = 0. 

Holevos uniqueness proof suffers from a (presumably only technical) gap. If, against expec- 
tation, it cannot be closed, let me point out that, at least under the additional assumptions 
in the next paragraph, the uniqueness of the strong solution of (2.8) is ascertained by its 
martingale property formulated in theorem 2.3 (if there were two strong solutions and if)? 
with the same initial value, \\tpt — ipt\\ 2 was a martingale, so ||^o — V'oll 2 = would imply 
EH^ 1 — ipf\\ 2 = 0, so ipl = if) 2 a.s.). Of course, the question whether there are additional weak 
solutions remains. □ 



2.2 General linear SDE with constant coefficients 
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2.2. Regularitat. The weak solution ipt is obviously strong if it does not leave the 
domains of Lj and K - or of another operator C with even smaller domain (then the operators 
in the weak form of the equation can be rolled back according to (K*<f>,ij>t) = {^),Kipt) and 
the scalar products with test vectors <p be dropped). For a proof of regularity, it turns out to 
be appropriate to formulate conditions in terms of such a C rather than subjecting Lj and 
K to additional restrictions. "C-strong solution" of a conservative equation of the form (2.7) 
then denotes an adapted, a.s. D(C)-valued process with continuous paths that satisfies the 
equation and, moreover, the estimates E||^|| 2 < E||^|| 2 and sup 0<s<t E[|C^ || 2 < oo for all 
t > 0. 

Theorem 2.2. Consider a conservative equation of the form (2.7) with D(K) C D(Lj) - 
that is to say, the operators have to satisfy 

(2.20) W L M 2 ~ 2Re W ^> = 

3 

for all ip E D(K). Let C be a positive definite self-adjoint "reference operator" with the 
following properties: 

• D(C) C D{K) Pi D(K*). 

• There exists an orthonormal basis ex, e%, ■ ■ ■ of S) consisting of elements from D(C) 
such that Vn € N : £V ||L*e n || 2 < oo. 

• Let P n be the orthogonal projection onto fj n := (ei, . . . , e n ). Then there are constants 
a, (3 > such that 

(2.21) 2Re (Cx, CP n Kx) + ^ \\CP n L jX \\ 2 < a (\\Cx\\ 2 + ||x|| 2 + 0) 

3 

holds for all n E N and x G $j n - 

• sup ngN ||CP n ^|| < \\Cx\\ for all x G D(C). 

Finally, let ip be a ^So-measurable, a.s. D(C) -valued random variable with E[|^[|^ < oo (as 
usual, we write \\x\\c '■= \/\\x\\ 2 + HCxIPJ. Then there exists a unique C-strong solution of 
(2.7) which satisfies the above mentioned estimates for all t > 0. This solution also satisfies 

MOM 2 < e at (E||C^|| 2 + at(EM\ 2 + /3)) . 

Theorem 2.3. With tp t denoting the solution from theorem 2.2, \\ipt\\ 2 is a martingale, in 
particular E\\ipt\\ 2 =~ 1 ' ' 2 



2.3. Comment on the form of the theorems. The conditions in the theorems appear 
to be formulated ad hoc. The authors comment that they take inspirations from "quantum 
stochastic calculus" - however, without further explanation, the conditions still look like a 
list of possible problems which one might encounter when dealing with these equations, but 
the real essence of which has not been inquired yet. However, the existence of weak solutions 
is ensured for many concrete applications - e.g. if fj = L 2 (M. n ) and the operators K and Lj 
are "built out ouf" differential and multiplication operators, then "usually" D = D* = 
will be suitable domains. Finding a reference operator C seems to be a more touchy business, 
which is not solvable by a general algorithm. In [15] Mora and Rebolledo have found such an 
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Solution theory for linear SDE with unbounded operators 



operator i.a. for the stochastic Schrodinger equation 




tracing back to Diosi ([5]), namely an arbitrary natural power of N := x 
operator which is self- adjoint on 



2 — t^A — id, an 




Even if this result looks satisfying, one has to mention the articles [14, 1], in which the authors 
have tailored a special method for this equation (transformation to the equation of the non- 
self adjoint harmonic oscillator) and found an integral kernel which provides an explicit strong 
solution for arbitrary bounded L 2 initial conditions. So one might guess that our restrictions 
are due to the particular ansatz rather than the actual behaviour of the equations. 



CHAPTER 3 



Product formula 



1. Linear SDE in M m 

Define a Wiener process £t on a probability space (f2,2t, P) and consider the m-dimensional 
SDE 

, s dX t =AX t dt + BX t d£t 

(3 - 1} x =x 

with A, B £ ]j mxm anc [ x £ IR m . The deterministic Trotter formula and the above mentioned 
result of [8] suggest that its solution flow C S) t can be obtained from the equations 

(3.2) dX t = AX t dt 
and 

(3.3) dX t = BX t d£t, 
resp. their solution flows 

(3.4) A s , t = e^ A 
and 

(3.5) B Sjt = e 
via the product formula 



(£ t -£ s )B-^^B 2 



(3.6) C ot = lim (Uen A B 



n-1 
k=0 

Here and in the following, products of operators are to be understood in the order of their 
action on a vector, i.e. 

n 

Hk := H n ■ ■ ■ H\. 

k=i 

In case that single random variables Co,t can be approximated in such a way, in contrast 
to the deterministic case the additional question arises whether the right-hand side in (3.6) 
is also a good approximation of (Co,t)t>o when considered as a stochastic process, i.e. one 
should also compare their autocorrelations or inquire the convergence of whole paths. For this 
purpose, it turned out to be useful to restrict to finite time intervals [0,T] and to replace the 
process (11^=0 &^ A Bk t k+i t )t>o, which makes the same number of approximation steps for all 

n 3 n 

times, either by a process (fn,T(t))te[o.T]i which is defined as f nj T (^T) := Ilfco e™ A Bi T j+i T 
on lattice points t S {0, T} and constantly extended in between, that is to say 

13 



14 



Product formula 



fn,T (^T + t) := f n , T (*T) for t = \T + t with r £ (0, £), or by a process <?„ iT which 
is defined in the same way on the lattice points and interpolated continuously according to 
9n,T (n^ + T ) := e rA Bk T k T+T Hf-o e™ A B j T j+i T . In order to avoid long case differentia- 

n 1 n J n 5 n 

tions, one can write this definition by use of the GauB bracket as 

LfJ-i 

(3.7) gn,T(t) := e( t ~ l ^ l ^ A B:nt:T t TT 6^%^ 



fe=0 



As will be shown, the paths of these processes converge to the ones of Co,, uniformly in the 
quadratic mean. Thus, or notion of convergence is slightly improved as compared with [8], 
where uniformity is missing. Moreover, we have noted that, on the analogy of Chernoff's 
theorem and also beyond the results of [8], the factors F(s,t) := e^~ s ^ A B s ^ appearing in 
(3.6) can be altered, as long as they keep solving (3.1) "for infinitesimal time intervals" when 
regarded as stochastic processes in the time parameter t with fixed s: 

Theorem 3.1. Let F : {(x,y) £ R 2 : y > x > 0} -+ L 2 (Sl,IR m) 
following properties: 

(i) F(s, t) = l + (t-s)A+ (& - £ a )B + o(t - s) or 

(i) ' F(s, t) = l + f* AF(s, r)dr + £ BF(s, r)d£ T + o(t - s) fur t -> a 

(ii) For A > and u > s + A, F(s, s + A) and F(u, u + A) are i.i.d. 



be a mapping with the 



Define / n>T (t) := uBo * F &, *±lT) fort £ [0,T], i.e. f n ?(t) 

I nt I i 

te [gT,*±iT[, andg r , T := F ([§ j J, t) H^J, *±T). 

a) (i) and (ii) imply Esup 4g [ T ] \\f n ,T(t) - C ,i || 2 "^°°> 0. 



Yl%lF{iT^T) if 



3+1' 



b) (i)' and (ii) imply Esup 4g[0T ] \\g n ,T(t) - C 0tt \\ 2 n ~ > °°> 0. 

Remark 1. (i) and (i') are not special cases of one another because f* AF(s, r)dr = (t — 
s)A + o(t — s), but in general only f* BF(s, r)d£ T = — £ S )B + o(y/t — s) holds. 

Remark 2. Under the additional assumption that Esup tg j Tj ||o(i)|| 2 — > for n -+ oo, 

the convergence of continues to hold true under condition (i)' and g n> T also converges 
under (i): Assume that the convergence formulated in the theorem and, in particular, the 
boundedness in the form sup ngN E sup <t<T \\fn,T 
proven. Then, since 



E sup \\g n .T{t) - /„,T(t)|| 2 =Esup 
te[o,T] t 



F 



<Esup|| F 



nt 
~T 

nt 



2 < oo and the same for g n ^ are already 
T 



,t 



n 

T 

n 



,t 



1 



one would only have to show lim r , 
e.g. for case (i) 



Esup||F 



f 

T 2 



nt 
T 



n 



II 



Esup||(t 



T 



<3^Pf + 3Esup(6-e 



t 

LfJ ^) 2 ||i?|| 2 + (l), 



nt 



1 ) U,t( 

2 Esup||/ n , T (t)|| 2 , 

t 

l) || 2 = which, considering 

2\„o 



n 



)A + (it-i\^\T_)B + o{t 



nt 
T 



n 



3.1 Linear SDE in M ; 



15 



boils down to Esup t (£ t — ^ntj z) 2 — > and is shown in the following proof from (3.11) on. 

Beweis von Satz 3.1. In order to check the plausibility of the theorem and as a techni- 
cal tool we first show that f Hj T and g n ^x are bounded in L 2 (S1, M mxm ) on [0, T] by constants 
not depending on n. A crude application of the triangle inequality for the L 2 (Q, M mxm )-norm 
of the single factors along with (ii) leads to 

Lf J-i \ 5 Lf J-i - - i 



E||/ n>T (i)|| 2 < e [J ll^r,*±lr)|| 2 = [J ( E ll F (^ r '^ r )H 2 ) 

\ fc=0 ' ' / fc=0 ^ ' 

I Tit j 



I A|| + \ / ''l +f n J 

and, since ^1 + -^j « e Cv/ ™, does not ensure boundedness, but this can be improved by 

reading the norm on W nxm as the Frobenius norm generated by (A, B) := tr(A t B) (other 
possible norms would be equivalent anyway) and replacing the triangle inequality by an exact 
computation. In case a) this leads to 

(3.8) 

LfJ-i LfJ-i 
E||/„, r (t)|| 2 <E J] \\F{ k T k -^T)f = nn k T k —T)\t 

Lf J-i 



n n ^ n n 

fc=o fc=o 



H (ni + -^ii 2 +E||(eMi T -efe r )5|i 2 +iE||o fc (-) n 2 + 

fc=0 v v 7 

+2E(1 + —A, (^fe+i T - £k T )B) + 2E((efc+i T - £fc T ) B > ( " J ) + 2E (! + ~A o fc { - 

77 n n ran y fi y 71 

Lf J-i 

^ IT ff 1 + T ||A||V + E(^ T -^ T ) 2 ||i?|| 2 + of^V 2E ^T-^T)( 1 + T A5) 

+ 2 (E|« ¥T - «j T )B|") 1 (E|, (i) f ) * + 2||1 + ^|| (e||, (i) II 2 ) ') 



1 + -IIAlA 2 + ^ IIBII 2 + o C -i +0 + 2./-||B||o (-) + 2111 + ^ Alio f i 



n" "/ n" \ if 1 ) V n" \n / n " \n 



Lf J-i 



< 1 + 



This term remains bounded because (1 + £) n — >■ e c . Now we have avoided a -4= summand 
within the bracket because E(l, {ik±\ T —ik T )B) = was computed instead of being estimated 

n n 

by Cauchy-Schwarz. The same trick works in case b), (£fc+i T — £k T )B now being replaced 

n n 

by J* ™ BF (^T, r) d£ T ; as for the other terms, it suffices to estimate the integrals via 
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Product formula 



J* a T dr = (t — s)a s + o(t — s) and f^a T d£ T = — £ s )ao + o{y/t — s) (valid in L 2 ) and 
products of them via Cauchy-Schwarz. 



The observation that / n ,T and g n ^r deviate only a little from the flow Co,* can be formulated 
in terms of integral equations similar to 



C ,t = 1 + f AC^ s ds + /* SC , s d6 : 
Jo Jo 



In Fall a), one gets 
(3.9) 



Lf J-i 



Lf J-i 



i+ E I? 



k=0 
If J-l 



n n 



r - 1 /n,T 



n 



,+ £ (^ + ({ ¥r -{ |T ) B + o l (I))/,„ r Qr) 



fc=0 

L^J-i / ,fc±i 



»+ E 

fc=0 

I rat I T 
L T J n 



Af n ,T [-T)ds + 



n 



k + l rp 



Bf n ,T 



-T J d£ s + o fe 
n / 



/n,T 



n 



Lf J-i 



= 1 + ^ ^/n,r(s)ds + ^ B/„,r(*)d& + E 0fe ("J ^ (^ T J ■ 

Using this representation, we aim at comparing f n T and Co,, by means of the Gronwall 
inequality. For this end, consider d n: T(S) := Esup 0<t<5 \\f n ,T(t) — Co,*|| 2 for arbitrary S <T. 
In 



d n ,T(S) <5E sup 
0<t<S 



+ 5E sup 

0<t<S 



I rat I T 
L T J n 



+ 5E sup 

0<t<S 



A(fnA s ) ~ Co,s)ds 



AC 0>s ds 

I rat I T 
L T J ra 

Lf J-i 



+ 5E sup 

0<*<S 



I rat I T 
L T J ra 



B(/n,T(s) " C 0)i )d& 



2 




+ 5E sup 


f BC 0:S d£ s 


0<t<S 





k=0 



fn,'. 



n 



the first summand can be estimated according to 

■L 



E sup 

o<t<s 



nt I T 
T J ra 



A(fn,T(s) - C 0tS )ds 



I rat J T 

2 <E sup L^J T / T " \\A(f n , T (s) - C ,s)\\ 2 ds 
o<t<s J n Jo 

s 



<T\\A\\ 2 E [ \\f n , T ( s )-C , s \\ 2 ds = T\\A\\ 2 f E\\f niT (s)-C , s \\ 2 ds<T\\A\\ 2 [ d n , T (s)ds, 
Jo Jo Jo 
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for the ensuing stochastic integral one additionally uses Doob's inequality and Ito's formula 
and gets 



E sup 

o<t<s 



I nt I T 
L T J • 



B(f n , T (s) ~ C 0)S )d6 



2 < 4E 



|"S|T 
L T J ; 



B(/n,r(s) - Co,,)d& 



- 4 / 

JO 



L t J " 



E||B(/ n , T (s) - C ,sWds < 4\\B\\ 2 / d n , T (s)ds. 

Jo 



In the remaining summands d n ^r does not appear - if they converged to uniformly in S for 
n — > oo, then 

d n , T (S) < (T\\A\\ 2 + A\\B\\ 2 ) / d n>T (s)d S + o(l) 

Jo 

would imply at first d n>T {T) < (l) e ( T !l^!l 2 + 4 ll B ll 2 ) T and for n — > oo the assertion of the 
theorem. 

Whereas 



E sup 

o<t<s 



\nt\T_ 
L T J n 



AC , s ds 



T rT 

2 <-|U|| 2 



n 



f nco,. 

Jo 



\ 2 ds = O 



n 



is straightforward and also the "little o" term, given that o& (^) is independent of f Uj T (^) 
and E||ofc(^)|| 2 = o(^) holds, satisfies 



Lf J-i 



(3.10) 



E sup 

0<t<S 

n-1 

<nJ^E 

k=0 

<o(l) sup E||/ n , T (t)|| 2 = (l), 

te[o,T],nen 



k=0 

Ok (-) fn,T (-T 



Lf J-i 



< no 



v 7 fc=0 



/n,T ~T 
n 



the remaining proof of limn-^oo E sup 0<t< 5 f,*nt , t BCo :S d^ s 
effort. It holds that 



requires somewhat more 



(3.11) 



E sup 

o<t<s 



\nt\T_ 
L T J n 



£C , s d& 



< E sup 

0<t<T 



\ nt\T_ 
L T J n 



5C 0)S d6 



=E sup 

0<k<n-l,0<t<^ 



Mkrp +t — Mkrp 



with the L 2 -martingale Mt := Jq BCo iS d^ s . First let n be fixed. By the Doob inequality, 
all increments Mk T+t — Mk T with fixed k and t G [0, ^] can be bounded by moments of 

n n 

Mm-i T — Mk T and increments over disjoint time intervals are independent. If M t were replaced 



n n 
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Product formula 



by the (one-dimensional) Wiener process £t, they would have identical normal distributions 
and all this would give rise to an explicit calculation: 



E sup 

0<k<n-l,0<t<± 



/ f 1 - P f Vfc G {0, ... ,n- 1} : sup 

Jo V V o<t<i 



sup 

0<fe<n-l,0<« 1 



n " 



> s ds 



< s ds = 



sup |6 - Col 2 < s 



ds 



1- 1 



< 



1 - 1 A 



ds 





1- 3 


^1 - 1 A 


TJ 2 S 2 



sup |6| 2 > s 



ds 



ds < 



In the penultimate step, one has to notice that, for small s, Doob's inequality gives negative 
lower bounds for 1 — P ^sup 0<t< i |£t| 2 > s^j , which have to be improved by the trivial estimate 

< 1 - P (sup 0<<< i |&| 2 > s) < 1. Now 

3 



lim E sup 

n->oo o<fc<ra-l,0<t<i 



lim 

n— >co 



1 - 1 A 



n 2 s 2 



ds = 



can be shown via the dominated convergence theorem: For the pointwise convergence of the 

integrand to 0, it is sufficient to show lim a ;_ > . 00 (l — = 1 or, equivalently, 

lim^oo — ^- In (l — ^) = for all c > 0. This follows from de l'Hospital's theorem because 

= -f^— — > 0. Moreover, for s > 1 the integrand is monotonically decreasing in n, 



_d_ 

dx 



so that m(s) : = X[o,i]0)+X(i,oo)( s ) (l ~ 1 A |l - J> |) = X[o,i] ( s ) + J?X(i,oo)( s ) is an integrable 
major ant. 

Note that the lattice points \T play no role and 

E sup ||&-6|| 2 <3E sup (||6-£|^|t|| 2 + ||£ S -A T ||2 + ||* -C,n.|T|| 2 ) < 

1 1 n 1 1 n 

9E sup ||6-6H 2 ^0fiirn^oo 

— — ' — — n 

holds as well. 

The m 2 components of the general L 2 -martingale M t can be represented in the form M\ )3 = 
C(^Mi,j) t as time-transformed Brownian motions (M' l, i)t = f Q \(BCo :S )ij\ 2 ds denoting the 
quadratic variation. For arbitrary M > 0, consider the decomposition 



E sup \\M t - M s 

\t-8\<Z 
1 1 n 



E 



Xsup 0<r<T \\BCo, r \\ 2 <M 



)+e(. 



■ Xsup < r < T \\BC ,r\\ 2 >M 



On the set considered in the first summand, \ {M lJ )t — {M l,3 ) s \ < M\t — s\ holds, so the time 
transform shows that the increments of M\ ,d can be controlled by the ones of a Brownian 
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motion and one gets 

E (•••X S up < r < T ||BC ,r|| 2 <M) < I 

E 



SUP \i( M i, j)t -t( M ij),\ Xsu Po < r < T ||BCV|| 2 <M ^ 



hJ 



E snp ™> ; 
l*-l<^ 



for the second summand there is no such control, so we estimate it at the outset by E ^sup| t _ s | <T 
According to Doob's inequality, one knows at least 



E sup \\M t -M s \\ 2 <4E sup ||M t || 2 < 4E (||M r || 2 ) < oo 

s|<T- J \0<t<T J 

and since the sets {sup 0<r<T ||-BCo, r || 2 > M} descend to the empty set for M — > oo because 
of the continuous paths of Co t, 



follows. 



The proof of part a) is therewith complete; in case b), a similar telescopic sum ansatz as in 
(3.9) yields 



(3.12) 



Lf J-i 



9n,T(t) = 1 + {an,T (^-^T^j - g UjT (^-T ) ) + g n , T {t) - 9n,T 



Lf j-i 



nt 



n 



nt 



t) ~ 1 ffn,T 



Lf J-i / ,,fe+i T 



J + E 



fc=0 



' I nt I T 

Lrjn 



Ag 



„, T (s)ds + J k 



k , 



' I nt I T 
L T J n 







) 5n,T f 





Lf J-i 



= 1 + jf Ag n>T ( 8 )d8 + jf J B 5 „ iT (s)d^ + °k 5n,T + O 5n,T ( 



nt 
T 



and one can continue along the lines of the preceding part. □ 
Corollary 3.1. Formula (3.6) holds true, the limit being uniform in the quadratic mean. 
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Proof. F(s,t) := e^ s ' A B s> t satisfies part b) of the theorem: (ii) folllows from (3.4) and 
(3.5) from the ltd formula one knows 

e^ A B s , t = l + f e^ A AB SiT dr+ I e^ A B s>T Bd^ T 

J S J S 

= 1 + f Ae^ A B s>T dr + I e( T -^ A BB s>T d£ T 

J S J S 

= 1 + f Ae {T - s)A B s , T dr + / Be^- s)A B s , T d^ T + ! \e {T - s)A B - 5e <T - s)A )5 S)T d( T 

•/ S *J S " s 



with 



IE II — 



1 



(t-s) 



— [\e^ A B - Be^ A )B s T d£ T || 2 = 
E\\(e (T - s)A B - Be (T - s)A )B ST \\ 2 dT fur t 



s 



because 



— E\\(e^ A B - Be^ A )B at \\ 2 <(t-a) 
t — s 



Ms)A _ i Ms)A _ i 

-B-B 

t — s t — s 



% • || AB - BA\\ 2 -1 = 0. 



This implies the convergence of g n ,T- Since i?o,t is an L 2 -martingale, the part of the proof 
starting with (3.11) tells us that even Esupi t _„i<t ||-£>s,i|| 2 ra ~" >00 > 0, so, according to remark 
2, the convergence remains true without the prefactors of g n ^T in the form (3.6). □ 

Theorem 3.1 also enables us to split off only a part of the deterministic term: 

Corollary 3.2. Let A = A\ + A 2 . Then C ,t = limn-»r» (Ilfc-o en Al D kt k +i t ) where D s>t 

n 1 n 

denotes the solution flow of dXt = A 2 Xfdt + BXtd^f 
Proof. Since 

e (ts)A 1Dg t = 1 + T e^-^AxD^dr + f e^ A2 A 2 D s , T dT + I e^ A BD s>T d^ T 

*) s J s J s 

=1 + f Ae^ A D s>T dr + I Be^ A D S:T d^ T + f\e^ A B - Be^ A )D s>T d^ T , 
the proof works as in corollary 3.1. □ 

Part a) of the theorem obviously comprises the Euler-Mayurama formula, the notion of con- 
vergence being slightly improved w.r.t. [13]: 

Corollary 3.3. Ul=o fl + ^ A + (£m t - £ it )s) ^> Cb )t . 

\ n n / 

Moreover, one can approximate the factors in corollary 3.1 by their first-order approximation: 
Corollary 3.4. UlZo (l + ±A) (l + (^ - £ At )fl) ^> C 0)J 



3.2 Equations with bounded coefficients 
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Proof. F(s,t) := (l+(t-s)A)(l+(^ t -^)B) = i+( t - s )A+(&-&)B+(t-a)(&-&)AB 
satisfies part a) of the theorem. □ 

As already the one-dimensional example dXt = (1 + l)Xfd£t shows, the stochastic part can in 
general not be further decomposed: One has Xt = Xoe 2 ^~ 2t , but the solution of dYt = Ytd^t 
w.r.t. the same initial condition is Yt = Xoe^*~^, so the solution flow is B S) t = e^' - ^ 2~ 
and rifc=d ( B k t+iBk k+i) X = e^^Xo ± X x . 



2. Equations with bounded coefficients 

The presented proofs remain true almost without change if one reads (3.1) as an equation 
for a process with values in a Banach space E, A and B as bounded operators on E, || • || as 
operator norm and the integrals as strong Bochner integrals or, respectively, their stochastic 
analogue. Only two parts have to be revised: 

1. The operator norm is no longer equivalent with a norm generated by a scalar product, as it 
was used in (3.8), so the application of the triangle inequality can no longer be circumvented. 
But there is another way to improve the estimate (1): Also by using the submulitplicativity of 
the norm, one destroys important information about the dynamics of the process - one takes 
into account the worst case that the oscillation in a time interval [|T, ^pT] amplifies the 
preceding one, although they typically cancel. One can start from the representation (3.9) 
and do already this step via the Gronwall inequality, which maybe does a better job in taking 
into account the described behaviour: Along the lines of the last section, one gets 



su P E||/ niT (t)|| 2 <4 + 4supE|| / T " A/ niT ( S )d S || 2 + 4supE|| / T " Bf n ^ T {s)A^ s f 

neN neN JO neN JQ 



Lf J-i 



+ 4supE|| V o k I - ] f n , T [ -T 



1\ . fk 



<4 + 4(T||A|| 2 + \\B\\ 2 ) f supE||/ n , T ( S )|| 2 d S + o(l) supEH/^)!! 2 ; 

JO neN neN 

in the last summand, a weakening of (3.10) was used. If N is chosen so large that o(l) < \ 
for n > N, then 

\ sup E||/„, T (t)|| 2 < 4 + 4(TP|| 2 + ||i?|| 2 ) f sup E||/ n , T ( S )|| 2 d S 

1 n>N JO n>N 

follows, so sup t6 [ T ] n>A rE||/ ra) j'(t)|| 2 < 8e 8T ( T ^ A \\ 2+ ^ B ^ 2 \ which is sufficient since we are in- 
terested in n — > oo. 

t 2 

2. The sum appearing in the proof of lim n _s. 00 

E sup < t < T /, nt I t BC 0tS dQ\ = would 
now be infinite and not obviously convergent. The claim can still be reduced to its one- 
dimensional form, namely via the integrator For adapted, E- valued processes of the form 
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at(uj) = J2i=o cl 2lt( uj )^( i-T i±iTl (*) anc ^ bounded by C > 0, for n > fc (that is to say, if no 

fc V fc 1 fc J 

more than one of the points xT lies between t and £ + ?) either 



a r d£ r 

(if no one lies in between) or 



a,M,T(6-6) <C"|&-& 

L T J fc 





2 




2 




2 


/ a r d^ r 


< 2 


/ a r d£ r 

J s 


+ 2 


/ a r d£ r 

■'I T 7 ! fc 


< 



< 2C [(£rfcjn t - £ s ) + (6 - ^rM-iT, 

I T I fc I T I fc 



<4C 2 sup {it- is? 

I*-I<£ 

(if one lies in between), in any case 



lim E sup 

| f _ s |<Z 
1 1 — n 



a r d£ r 



< lim 4C 2 E sup (6 - 6 







1 1 — n 

holds. As such processes at are dense in L 2 (Q,E), it only remains to show the continuity of 
the limit in at- This can be done via Doob's inequality:: 

2 

lim E sup 



1 1 77. 



/ (a r - a r )d£, 

J s 

< lim E sup 2 (a r — d r )d£ 

\ts\<^ V JO 

/ (a r — d r )d£ r 

JO 



+ 



<4 lim E sup 

n ^°° 0<t<T 



< 16E 



J (a r - a r )d£ r \ 
/ (a r — a r )d£ 



<16E 



|o r — d r \\ dr 



3. Equations with unbounded coefficients 

As outlined in chapter 2, the stochastic part of an equation like (2.7) with unbounded op- 
erators should not be considered separately, but rather in connection with a deterministic 
correction term in the form (2.3). Therefore a product formula analogous to (3.6) can be 
most easily found for the equation (2.8), which explicitely contains this corrector term. 

The proof of the result formulated below shows in particular the weak solvability of (3.13) 
via the product ansatz (3.14), which is at the same time a more explicit representation of the 
solution than in [10]. We use ideas from [10], but our proof is independent. Also we have 
to formulate ad hoc conditions, but they appear to be rather natural: Roughly speaking, we 
require the existence of an invariant subspace which the partial solution flows Ht and A Sj t do 
not drive out of the domain of each other's generator. Such conditions appear in deterministic 
semigroup theory as well, the only difference being that a well-established terminology permits 
more handsome formulations. 



3.3 Equations with unbounded coefficients 
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Theorem 3.2. Consider the stochastic Schrddinger equation 

(3.13) dip t = (-iH - ^A 2 )^ t dt + Afodtt 

with self-adjoint operators H and A and the partial solution flows H t := e ttH and A Sit := 
e (£. t -£.s)A-(t-s)A 2 _ Suppose that there exists a dense subset M. C S) such that for ip G M and 
sufficiently large n all products appearing in (3.14) (also the empty one) are in Dom(A 2 ), 
A t t+ T \\ . . . € Dom(H) for all t and Ht \\ . . . € Dom(A 2 ), and another dense set M C f) 

such tha for all $ e N both lim^o AHf(f> = A(f> and lim^o A 2 Ht<fi = A 2 (j) holds. Then the 
equation is weakly solvable and, for all T > and ip £ f), its solution flow Co,t is given by 

Lfj-i 

(3.14) Co ti> = lim n HTAk T k+i T tfj 

' n-^OO J- J- n n > n 

k=0 

with weak limit in L 2 (f2,f>). 

//, in addition, the conditions from theorem 2.2 about regular solvability are satisfied and 
M.C\D(C) ( C denoting the reference operator used therein) is dense in $), then the convergence 
is even strong. 



Proof. Choose ip 6 M and let ip be a product appearing in (3.14). The conditions make 
sure that, for s > 0,t € (s,s + the Ito formula can be applied to H t _ s A s jip and, doing 
this, one gets 

Ht- a A, tt j> =Y> - ij H T _ s HA S:T 4>dT + J H T _ S A S ^A^ T ~\ J H T ^ s A SjT A 2 ^dr 



-iH - ^H T - s A 2 H s _ T j H T _ s A StT 4>dT + jf {H T _ s AH s - T )H T - s A s , T i>d£ T . 

I n* 1-1 

With the ansatz pursued in (3.12) for g n ,T{t) ■= H.\nt\TA\nt\T t Yl^Q HrBk T k+i T 

L T J n L T J n ' n n'n 

one gets 

Lfj-i 



n n n \H J 



(3.15) 



k=0 

+ {^H t _ ynt\T_A \ nt | T 



Tin LTjn 





nt 


) 9n,T ( 


J 



n 



,t u 



= -i Hg n)T (s)ds+ I # s _i ™ | ™ | z_ s 5n,r(s)d£ s 

JO JO IT J n L T J n 

1 ^ 

^- L¥ j^ 2 ^L¥JF^ (s)ds - 
Ths is a linear equation for g n ^ with self-adjoint coefficients 

J^(n>,t) ■= H t _ynt^TAHyrit^T_ t Und ^4 2 := ^ntj Tj4. 2 -ff|^ratJ 
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which, since < s — |_^rj ^ < ^, converge uniformly in t to their counterparts in (3.13) on 
Af for n — > oo. Thus, we are in a similar situation as in the proof of theorem 2.1 (after having 
arrived at (2.12)) and, like there, can pass over to the weak limit: 

Since H t and A st preserve the L 2 (n,fi) -norm, 

(3.16) E\\g n , T (t)i;\\ 2 = U\\ 2 

holds, so the terms g n ,T(-)ip, viewed as mappings from a time interval [0,T] to L 2 (Sl,fj), 
are uniformly bounded for any ijj G ft. We would like to show equicontinuity and apply the 
Arzela-Ascoli theorem. Since we have no clue concerning the limiting behaviour of expressions 
like KWA^'^ g n) T(t)^>\\ 2 , we pass over to weak convergence in L 2 (f2,f>) = L 2 (f2,C) (8>.f), i.e. 
we choose ^ 6 L 2 (Q,C) and € AA and consider the complex- valued functions f n (t) := 
^(<f>-> g n ,T(t)ip}- Due to the boundedness of g n: T{t)ip, we only have to consider £ and from 
dense subsets of the unit spheres, therefore we restrict to <f> € M; since L 2 (f2,C) and f) 
are separable (the first statement follows from the Wiener chaos decomposition), even an 
appropriate sequence (6^)ieN suffices. The f n (t) are uniformly bounded according to the 
Cauchy-Schwarz inequality and (3.16) and equicontinuous according to 

int)-n*)i 2 <m 2 -m,9n,Tm-9n,T(s)^)\ 2 

=E| £<0, A^g n , T (r)^ T - j\cf>, (iH + l -A 2M )g n ^)dr\ 2 

<2E| /"*(i (n ' T) *0, ffn,T(r)V)d£ T | 2 + 2E| f {(iH + i 2 ^*)^, ^(r^dr) 2 

is J s 

<2 /"*E| (1(^)0, ff „, T (r)V)| 2 dr + 2|t - s\ [' E\{(iH - i^V, <? n , T (r)^)| 2 dr 
<2|*-a| sup ||I("' T V|| 2 sup E|| 5niT (r)^|| 2 

re[s,t] re[s,t] 

+ 2\t-s\ 2 sup \\{iH - A 2in,T) )4>\\ 2 sup E|| 5riiT (r)V|| 2 

re[s,t] re[s,t] 

<2\t-s\ sup ||I( n ' T V|| 2 + 2|t-s| 2 sup \\{iH - A 2{n,T) )<j)\\ 2 

re[s,t] re[s,t] 

<c^ T \t-s\\m 2 . 

Therefore, there exists a subsequence (nj^keN of (n) ne jq such that E^-^i, g n i T (-)-0) converges 
uniformly in t, a subsequence (n|) of (n£) such that also E£ 2 (02j <7 n 2 t(')V') converges and 
so on. For the diagonal sequence (n^) =: (n&) obviously E£j((j)j, g n k T (-)ip) converges for all 
j £ N uniformly in t; since 

L 2 (n, S))3r]^ lim E(r?, g n k T (t)i/j) G M 

fc— S>00 fc' 

is a continuous functional, there exists a process V't with 

lim E£(0, 5n j iT (t)^) = E£<0,^> 

for all £, </> uniformly in t. 



3.3 Equations with unbounded coefficients 
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In order to show that this is the wanted solution, the convergence of the remaining summands 
in the "weak form" of (3.15), i.e. in 



(3.17) 



m<t>,9n k , T m) =M(^) + E£<0, / A^g nktT {s)^ s ) 

Jo 

- E£<0, J (iH - ^ ink ' s) g nktT (s)^ds) resp. 

m<f>,9n k ,TM) =m^i>)+^ Ai K ' S) 0,5n fc ,T(s)^>d6 

J 

+ El j\(iH - l -AZ {nk ' S) )<P ignk , T {s)^)&s 



(strong Bochner integrals and their stochastic analogue allow for dragging bounded operators, 
in particular scalar products, inside), to their counterparts in the weak version of (3.13) has 
to be checked. To begin with, 



mj ((iH- ii 2( " k ' S V,9„ t ,r(#d S -E^ {{iH - l -A 2 )^ s )ds\ 

=\J*m(iH- ±A* (nk ' s) )<i>,g nk , T (sW)ds- J*^((iH- ±A 2 )<t>,i/> 8 )ds\ 

<\ j\l(\(A 2{nk ' S) cf> - A 2 <P),g nk , T (s)4;}d S \ + | f^l((iH - \a 2 )4>, (g nk , T (^ ~ ^s))ds\ 

<(E|£| 2 )§ sup ||i(i2 (nfc ' s V-A 2 0)|| f (E\\g nktT (s)^\\ 2 ^ds + ...^0; 
se[o,t] z J o 

for the summand in the middle of (3.17) we restrict £ WLOG to the iterated stochastic 
integrals Ii{t) specified in lemma 2.1. This restriction to a total subset is justified by the fact 
that, according to 

E| /V" fc ' S V,2n fc ,TW>d6| 2 = fn(A^<P,g nk , T {sm 2 ds 
Jo Jo 

< sup \\A^ s Uf fng nk ,T(sWte<CT 
se[o,T] Jo 



the sequences of stochastic integrals which are checked for convergence are bounded. Then, 
similar as above and in addition using the Ito formula, the independence of the Wiener 
processes and the fact that, according to their recursive definition, all Ii(t) are L 2 -martingales 
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and have mean value 0, we conclude 

Jo Jo 



=\E /*a( S )/,( S )d^ s f\A^ s U,9n k ,T(s)^ s -E/*a(s)/Ks)d6 f (Acf>^ s )^ s \ 
Jo Jo Jo Jo 

=\E f^(s)W)(^ nk ' S U,9n k ,T(^)ds-E f a{s)Ii(s)(A(/>,^ a )ds\ 
Jo Jo 

=| f Ea( S )/ i (s)((I (n '=' s V - M g nk , T (s)i;}ds\ + | /' Ea(s)J,(a)(A& (g nk , T (^ ~ ^s))ds\ 
Jo Jo 

<(E|/ ; (s)| 2 )^ sup \\A^> S U ~ A0W f E\\g nkiT (s)i;\\ 2 ds + . . . ^ 0. 
se[o,T] Jo 

Therefore, the constructed process tpt solves equation (3.13) (and agrees with the solution 
from 2.1 insofar as the solution is unique). The same method also allows to start with any 
subsequence gn k ,T(t)ip and extract a sub-subsequence converging in the same manner weakly 
to ipt, so even g n ,T(t)ip — 1 ipt holds true. 

In* |_l 

Compared to f n ,T(t)tp '■= Ylk=o HrA kT k+i T ip, the product representation (3.14), g n ^< 



n n ' n 



contains some additional prefactors, but they do not matter because 

Lfj-i 

PW, fn, T m - 9n, T m)\ 2 = \^{4>, {H t _, M l T A { ^ l T t - 1) 17 HtB^ k+i^ 

L T A n L T J n 1 n n ' n 

k=0 

Lfj-i 

= |^<(^|nl|T f H* -1)0, IT HrB kTk±lT ij)\ 2 

LtJu' 1 1 LtJt. - LJ - ii n ' n 

fe=0 

Lfj-i 

<E|e| 2 E||(^ LfJfit ^ +LfJ? -l)0|| 2 E|| n 



fc=0 



=E|e| 2 E||(^ Lf jz )t ff_ t+L ^ j z - 1)0H 2 E||V|| 2 -> 0, 

since L^J ^ — >■ i and A s j and -ffj, regarded separately, are strongly continuous and uniformly 
bounded because of E||A Sit VI| 2 = ^ll-f^ll 2 = IIV'II 2 (for A s>t see section 2.1). 

Since lin^-^ E||/ ni T(t)V'|| 2 = HV'II 2 holds true and, under the conditions of theorem 2.3, 
ip G M Pi D{C) also satisfies E||-(/>t|| 2 = HV'II 2 , m this case the convergence / n ,r(i) — > Co,* 
is even strong. Since Co,t as well as all operators f n ^r have L 2 (0,,£))-norm 1, the strong 
convergence on the dense subset M n -D(C) implies the convergence on the whole of ft. □ 



Certainly it would be desirable to prove uniform convergence in the quadrativ mean as in 
theorem 3.1, however, the way to get there is less obvious in this case than in the ansatz 
using the Gronwall inequality. 



3.3 Equations with unbounded coefficients 
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In order to get an idea how one can find M. and M in concrete cases, let us again have a look 
at Diosi's equation 

(3.18) dip t = ~\Mt + s^d& - ^arVtd*. 

We try the maximally possible choice 

M = AT = D(A) n D(x 2 ) = {0 G H 2 (R) | - 2 0(-) G L 2 } =: D. 
U c/) eT>, then obviously Aj,t0 G f, since 

and multicplication with a Gaussian preserves decay and smoothness properties. H t (j) G D(A) 
ist well-known and instead of H t <f) G D(x 2 ) it is easier to show H t 4> G _ff 2 (]R) (Sobolev space): 
This follows from the fact that H t <fi solves the Fourier transformed Schrodinger equation 

i/jt(0 = ~ |s an d is therefore given by H t 4> = err *0. It remains to be shown that 

T> is also a suitable choice for M\ according to the Plancherel formula, one has to check 

xH t (j)(x) x(j)(x) and x 2 H t (/)(x) x 2 (j)(x) for t — > or, equivalently, ^H t <p(^,) J:<A(£) 
and ^Hrf>(0 Jp<K6- The formula for H^> tells us 

= e-^ H {-itm) + |<K6) and 

^ w«) = e-^ H \-{t 2 e + it)to - 2«^(0 + ^toi 

and the assertions follow from the dominated convergence theorem since <f> G H 2 and ■(/>(■) G 
i/ 1 imply that 2 4> and ■^<^(-) are in L 2 . 

The fact that T> agrees with the domain (2.22) of the mentionned reference operator implies 
strong convergence in the product formula. 

The application in chapter 4 requires the following straightforward generalization of theorem 
3.2: 

Corollary 3.5. Consider the stochastic Schrodinger equation 

m ^ 

(3.19) dipt = -iHip t dt + 2^(A^ t d# - -A 2 ip t dt) 

with independent Wiener processes self-adjoint operator H and commuting self-adjoint 
operators Aj and the partial solution flows H t := e ltH and 

j 

Suppose that there exists a dense subset A4 C ij such that for ip G M. and sufficiently large 
n all products appearing in (3.20) (also the empty one) are in P) Dom(A 2 ), A..,t \\ . . . G 

Dom(H) for all t and Ht Y\ ■ ■ ■ G f] Dom(A 2 ), and another dense set M C S)such that for 
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all (j) 6 M and j both lim^o AjH t (j) = Aj(f> and lim t _^o A 2 H t (j) = A 2 cp holds. Then the equation 
is weakly solvable and, for allT > and tp £ $), its solution flow Cb,t is given by 

Lfj-i 

(3.20) C ot ip=\im TT HrAkrp k+irpip 

n— >00 J-- 1 n n ' n 

k=0 



with weak limit in L 2 (f2,f)). 



//, in addition, the conditions of theorem 2.2 about regular solvability are satisfied and Ai n 
D (C) (C denoting the reference operator used therein) is dense in S), then the convergence is 
even strong. 

Proof. As in the proof of theorem 3.2 one gets by means of the ltd formula 

H t - s A Sit ij) = 

=^- i j a Hr-sHAs^dT + ^2(f s ffr-A,T^^? ~\j 8 H r-sA s , T A 2 4>dT) 

=4> + J -iHH T _ s A s ^dr 

+ (H T - a AH 8 - T )H T - 8 A aiT $dZ T - X - j (H T - s A 2 H s _ T )H T _ s A S:T 4>dT. 

j J S J S 

(Here it is important that the Aj commute - otherwise, terms like 

k— 1 m 

TJ e (£-m-(t- S )A 3 2 ^ "Q e (£i-ti)A j -(t-s)A* 
j=l j=k 

would appear and would not agree with AkA s> t because A^ could not be put in front. In order 
to deal with such cases, one would not only have to require compatibility conditions between 
H and Aj, but likewise for the Aj among themselves.) 

Via the same telescopic sum ansatz one gets the equation 

9n,T(t)ip = - i J Hg ntT (s)ds + ^2(J H s _^^AjH^^_ s gn,T(s)d£ s 

H s-W\l A2 3 H W\l-s9n,T{s)ds), 

the limit of which can be performed as in the previous proof. The £ needed for testing weak 
convergence are still provided by lemma 2.1. □ 



Equation (4.5), the multi-dimensional version of Diosi's equation (3.18), is covered by this 
corollary in the same way as (3.18) by 3.2. 



CHAPTER 4 



Example: Quantum mechanical collapse models - GRW and 

QMUPL 



In the end, I would like to present the two quantum mechanical collapse models which made 
me think about the product formula. Both models complement the conventional Schrodinger 
evolution fa = e~h tH (fi by a - mathematically well-defined - mechanism leading to a random 
spatial localization of the wavefunction und thus solving the "measurement problem" which 
is often illustrated by "Schrodinger's cat", fa is replaced by a stochastic process (we will 
not change the notation, though) on a probability space (fl,^, P). Arbitrary Schrodinger 
operators H are admitted, though, for simplicity, we will restrict to the case of a single free 
particle and omit purely physical constants, so ip G L 2 (M 3 , C) and H = — ^A. 

1.1. GRW. In the Ghirardi- Rimini- Weber model (shortly GRW, see [6], [2] and, for 
a mathematically rigorous treatment, [19]), the Schrodinger evolution is interrupted at the 
random jump times t\,t%... of a Poisson process with intensity \i (so, in paricular, Er^ = ^) 
by a collapse mechanism - namely, the multiplication by a Gaussian "hitting function" with 
spread For mathematical simplicity and - see later - in accordance with our needs - we 

replace the by deterministic times ^ Then the mechanism can be described as follows: 
We start with an initial wavefunction fa £ L 2 (IR 3 ,C), construct a process ((f>k) of wavefunc- 



tions on the space Q := (IR 3 )'-'"^''" ' , the coordinate projections of which will be denoted 
Y\ , Y2 , . . . , recursively via 



1. The models and the question 



1 2 



(4.1) 



ip n (Yl, ■ ■ ■ ,Y n ,x) : 




ll^n||x 2 (K 3 ) 



and define 



(4.2) 
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which reads explicitely as 



(4.3) 



\AY 1 eA 1 )= ! \\Myu-)\\ 2 L2 m ^yi 

J A\ 



J e- a ^ 2 \(e-i H <f> )(x)\ 2 dxd yi , 
J{Y 1 eA 1 } 

J{y 1 &a 1 }Ja 2 \^J \m\\ 



GRW 



« f e- a \ x -^ 2 \(e-i H ^)( yi ,x)\ 2 dxdy 2 - 1 \\Mvu ^fdyi 

a 1 ja 2 v j ImmrW 



a I _„,u_„„i2 
A 1 JA 2 * 



e «l a: -wr|( e -^ ff e -fl-Kil a e-i ff ^ )( a; )|2 da . dltedl/1 



and, inductively, 
(4.4) 



= / *P / |(e-fl-^l 2 e-^---e-fl-^l 2 e-^ ( /»o)(x)| 2 dxd(y 1 ,...,y n ) 

and yields a consistent family that has a projective limit on the whole of Q. The GRW process 
is obtained by canonically identifying Q, with a subset of C([0, oo),M 3 ), pushing P Qj/i forward 
to this space via the corresponding inclusion and extending (cp t ) according to <f>k, := e" 

between the collapse times. 



1.2. QMUPL. In the QMUPL (Quantum Mechanics with Universal Position Localiza- 
tion) model, tracing back to Diosi ([4], [5]), a time-continuous collapse of an initial wavefunc- 
tion 0o £ L 2 is constructed by the aid of the stochastic Schrodinger equation 

^ ^ dip t {x) = -iHifitdt + \f\ipt{x)x ■ d& - ^\x\ 2 ^ t {x)dt 

with a three-dimensional Wiener process (£t) on a filtered probability space (Q, (fit), Q) and 

an intensity parameter A > 0: The "physical" collapse process is 0t := jmr 1 > weighted 

HWIl 2 (r3) 

by a new measure defined by 

(4.6) F X (A) := E Q ( XA \\M 2 ) &t A G &. 

The martingale property of HV'ill 2 (theorem 2.3) ensures that this is indeed a consistent 
definition of a measure on f2. One physical justification of this model lies in the fact that the 
isolated collapse process (i.e. (4.5), (4.6) with H=0) is a continuous extension of the GRW 
collapse process (see the after next paragraph). 



4.2 Equivalence of the collapse processes 
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1.3. The question. In [4], Diosi has claimed that the QMUPL model can be obtained 
as a scaling limit from GRW by increasing the collapse frequency and, at the same time, the 
spread of the hitting function, thus weakening the effect of a single collapse and keeping the 
overall effect constant. The appropriate scaling of the GRW process for a — > was found to 
be 

(4.7) p = ™ 

a 

with constant A > 0. Diosis proof lacks mathematical rigor and only considers the level of 
the Lindblad equations for the statistical operator p t := E(\ip t ) {ipt |)> which are (in "position 
representation" ) 

d_ 
d~t r 

for GRW and 



- L Pt(x,y) = -^[H, Pt ](x,y)- fx(l-e y)2 ^j p t (x,y) 



—p t (x,y) = - % -[H,p t ](x,y) - ^(x - y) 2 p t (x,y) 
for QMUPL. For a — > 0, the GRW Lindblad term can be linearized in a and in fact one gets 

p (l - e-^ 2 ) « „ (l - (1 - SL {x - yf) =^ {X - yf. 

We have paved the way for proving the following result directly on the level of the wavefunc- 
tion: 

Theorem 4.1. For a — > 0, the finite dimensional distributions P Q 2a o (4>ti, ■ ■ ■ , (fan) °f ^ e 

' a. 

scaled GRW processes converge weakly to Pa°(<^d • • • j0t„) -1 ; which are those of the QMUPL 
process. 



It seems plausible that, apart form technical complications, random GRW collapse times 
would not change anything. Moreover, it is worth mentionning that the only purpose of 
the random collapse times is to avoid an arbitrary assignment and, in the case of several 
particles, an arbitrary coupling of the collapse times for the single particles. This arbirariness 
is removed by the continuum limit anyway. 



2. Equivalence of the collapse processes 



The comparison of GRW and QMUPL is particularly elucidating if one neglects the Schrodinger 
evolution (i.e. puts H=0). Then, according to (4.1) and (4.4), 

(4.8) </>n(x) = ce"f ^- yi l 2+ - + l x - y "l 2 ]0 o (x) 

(c is the Y-dependent normalization factor) with 
(4.9) 

F a ^(Y 1 eA 1 ,...,Y n eA n ) = x f^ [ [ e^-^) 2+ --- + ^) 2 ^o(x)| 2 dxd( yi ,...,y n ). 

The solution of the collapse part of (4.5) is 

^t(x) = e >/Ax*t-A|x| a t 
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and can be rewritten as 

ipn(x) =e f m 4>o{x) = e v m m mm <po(x) 

or, in normalized form and taking into account the scaling (4.7), 

(4.10) (f)n{x) = ce v ^ </>o W- 

The quantities Zf,, := -^=(^ — £ fc -i) formally correspond to the from GRW and, by the 

aid of (4.6), their distribution under ¥\ can be determined if one takes into account that they 
have independent centered normal distributions with Y(Zk) = 4A~ = 2a under Q: 

P A (Zi € A u . . . , Z n € A n ) = E Q (x{z 1 eA 1 ,...,z n eA n }\\^ f) 
= e q(X{z 1 6A 1 ,..,z„6A ?1 } / e " » \<po\ {x)dx) 

(411) =Eq(x { z i6 a 1> ...^ 6 ^} / ^^-^l^l^l^xjdx) 

= / / e EJ=iP«*- fc -a|x| a ]|^|2 (a . )da . Ke-aEJ^ l^l 2 d (^, . . . , z n ) 

= A /^ / /e-^^l^l^ol'^dxd^,...,^) 

(4.10) and (4.11) agree with (4.8) and (4.9). Thus, there is no qualitative difference between 
GRW and QMUPL processes with parameters (a, ^) and A - the first can simply be obtained 
by restricting the latter to the appropriate discrete instants of time. 



3. QMUPL as continuum limit of GRW 



The circumstance that the GRW collapse is a discretization of the QMUPL one suggests to 
read the iterative definition of the complete GRW model (including the Hamiltonian) as a 
product formula approximation of the solution of (4.5) - the only question to be settled is how 
the change of measure (4.6) fits into this framework. To this end, we not only approximate 
the solution of the QMUPL equation up to a time T, according to theorem 3.2, by 

'yj 1 >A(-)-K* ± i-*i)-£H a _l Hm 

J [ e MM g M (J)q ? 

k=0 

but also the change of measure by 

F X (A) :=E Q ( X A\m 2 )Kr Ae$ T . 
Then, withZ fc :=^=(£i-£^i), 



(4.12) 
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with ^-dependent normalization factor c and, as a generalization of (4.11), one finds that 
P£(Zi €A 1 ,...,Z n €A n ) = E Q ( X{ZieAl _ Zn€An} ||^ || 2 ) 

/n—1 

k=0 

= [ [ | JJ e a ^ Zk + 1 -^%~^ H (t)\ 2 (x)dxJ^ n e- a ^kZo l*fc+il a d(zi, ...,z n ) 

=\[^ I I lIIe^^-tl-l'-fl^^e-M^^^dxd^,...,^) 

/ — n n—1 

=\ - / / in e " fi ~ zfc+i|2e_ ^V(^d(^,...,z n ) 

V A" JAix---xyl n J fc=0 



if ^ < 7" ((•) denotes the argument of the function). This expression agrees with the GRW 
distribution (4.4), i.e. on 3r it holds that P^ = P Q)M with a = The question of convergence 
in theorem 4.1 for a — >■ resp. /i — > oo can now be settled without further reference to the 
GRW model by showing 



lhn E( X{0 M 6Ali ...^ 6An} ||V;fJ| L 2) =®(X{<t> tl eA 1 ,...,<p tn eA n }Ht n \\L*), 



tpt denoting the solution of the QMUPL equation (4.5). However, this is obvious because, 
since tp£ — > tp t 111 L 2 (Q, L 2 (M 3 )) according to corollary 3.5, the sequence of the integrands is 
convergent in probability and uniformly integrable. 



Note that the reverse order of the factors in the definition (4.12) of ipjf w.r.t. theorem 3.2 
plays no role: 



Akrp k + lrpHTtp 



I nt I -| I Tit I 

It J It J 

] f HrAk T k+± T 1p — I J 
k=0 k=0 
LfJ-2 

<^-i)^ LfJ _ 1)LfJ n h iHt,^ + 

k=0 




Akrp k + lrpH Tip 

n. ' n n, 



A q t(1-Ht)^, 
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Example: Quantum mechanical collapse models - GRW and QMUPL 



for the second summand, the unitarity of Ht and (2.5) imply E|| • || 2 = ||(1 — Ht 

n 

n — > oo; as for the first, we know at least 

[f J-2 

|Ef(0, (Ht - I) A I ut I , I M i TT HrAk T k+i T ^)\ 2 

n ITJ 'LtJ - LJ - n n ' n 

fc=0 

=|E« n a; Hi L4t u| _ |al| (ffi-i)^,«| 2 



\fe=LfJ-2 




<E|£| 2 E| 



n ^ 

v fe=[f J -2 

=E|e| 2 ||( J ff_T-l)0f 11 ' 112 



Ji k + l rpH T 

n ' n n 



(#_T - 1) 



0. 



Thus, despite the different order, ^ converges weakly and, since E||V>f || 2 
also strongly to V't- 



for 



Eli 



J t\\ , 



CHAPTER 5 



Conclusion 



We have found a new access to the solution of equation (1.1) by means of a representation 
which, compared to earlier works, is more explicit and closer to the physical motivation of 
such equations. Possibilities for technical improvements are obvious and werw mentioned on 
the spot. The applicability of our product formula to the QMUPL equation subject to the 
nonlinear modification by a measure change, as described in chapter 4, raises the general 
question to what extent such product formulas remain valid for diffusion processes described 
by nonlinear equations (in the deterministic theory of nonlinear semigroups, this topic is 
inquired e.g. in [12]). In particular, since the QMUPL model, as it seems, does not admit 
adaptions to further physical requirements (such as identical particles, Lorentz invariance, 
relativistic interaction etc.), an extension of the product formula to models like CSL ([7]), 
which are considered as more promising in this respect, is suggested. 
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